arXiv:astro-ph/0503472vl 22 Mar 2005 


EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 


CERN-PH-EP /2004-076 
December 8, 2004 


Measurement of the Shadowing of High-Energy Cosmic 
Rays by the Moon: A Search for TeV-Energy 

Antiprotons 


The L3 Collaboration 


Abstract 

The shadowing of high-energy cosmic rays by the Moon has been observed with a 
significance of 9.4 standard deviations with the L3+C muon spectrometer at CERN. 
A significant effect of the Earth magnetic field is observed. Since no event deficit on 
the east side of the Moon has been observed, an upper limit at 90 % conhdence level 
on the antiproton to proton ratio of 0.11 is obtained for primary energies around 
1 TeV. 
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1 Introduction 


Physics motivation The effect of the Moon, or the Sun, on cosmic rays was first noted by 
Clark in 1957 [1]. As these bodies pass overhead they block the particles, so their shadows in 
the cosmic ray flux must be visible by detectors on Earth. 

However the first observation of such a shadowing had to wait for the results of the CYGNUS 
collaboration in 1991 [2], There are two reasons for this long delay. First, the particles must be 
insensitive or weakly sensitive to the Earth magnetic field. Obvious candidates are 7 -rays or 
energetic cosmic ray particles. The former are very rare and the observation of the latter above 
the nearly-isotropic large background of low-energy cosmic rays, was only possible with the 
advent of large Extended-Air-Shower (EAS) detectors, able to collect high statistics. Second, a 
crucial parameter is the angular resolution of the detectors. The signal over background ratio 
is inversely proportional to the square of this angular resolution and events are spread out from 
the expected position due to the finite angular resolution. The performance of the detector has 
to cope with the angular radius of the Moon (or the Sun), each having approximately a 0.27° 
radius, and only at the beginning of the 90’s, the angular resolutions of cosmic ray detectors 
reached the one-degree level. 

Since then, several other experiments, both EAS arrays and large underground detectors have 
been able to see the Moon-shadow effect [3-7]. The observation is used for a check of the 
angular resolution of the apparatus and, by comparing the observed position of the deficit to 
the expected Moon position, to evaluate systematic pointing errors. The understanding of the 
alignment and of the angular resolution is a key issue for any point-source search. 

In 1990, a more challenging use of the Moon shadow effect has been proposed [ 8 ]. The use 
of the Moon collimation, together with the Earth magnetic field, allows a charge determination. 
Negatively charged primaries are deflected towards the west and positively charged primaries 
towards the east. If antiprotons are present in the cosmic ray flux, they will generate a shadow 
on the opposite side of the Moon relative to the shadow from cosmic rays induced from matter. 
This article discusses a search for cosmic ray antiprotons using the L3+C muon spectrometer 
of the L3 detector at the CERN LEP accelerator. Published direct measurements of the p/p 
ratio exist only below 40 GeV and this method is sensitive at TeV energies. Only non-standard 
sources would be the origin of such high-energy antiprotons. No such study has yet been 
published using the Earth-Moon system as a spectrometer: EAS arrays (with the exception 
of the Tibet array) and underground detectors have a too-high detection threshold, so that 
the effect of the Earth magnetic field is just a small perturbation. This is no longer the case 
in the L3+C experiment. Due to only 30 m of overburden, multiple scattering remains small 
even for low energies and the accumulation rate is much larger than in other underground 
detectors. Moreover, the measurement of the muon momentum allows for an off-line tuning 
of the threshold, leading to a possible optimisation of the shadow effect. The present study 
exploits these possibilities and looks at the conditions to set a limit to the p content in the TeV 
region. 

This section describes the status of cosmic ray antiproton data and summarizes other exper¬ 
imental observations of the Moon shadow by cosmic rays. The experimental setup is presented 
in Section 2. Section 3 contains a review of the main parameters involved in a Moon-shadow 
experiment and describes the role of the different Monte Carlo simulations. In particular, the 
simulation of the experimental angular-resolution is checked with the help of two-track events. 
Data and backgrounds are presented in Section 4. The observation and interpretation of the 
deficit of events in the Moon direction are described in Section 5. Results concerning the ex- 
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perimcntal angular-resolution and a possible shadowing effect due to antiprotons are discussed. 
Conclusions are presented in Section 6. 

Cosmic-ray antiprotons The experimental p/p ratio below 50 GeV is compatible with a 
secondary origin of the cosmic ray antiprotons [9]. Data are obtained by balloon-borne ex¬ 
periments and, recently, satellite experiments. The CAPRICE [10] and HEAT [11] collabora¬ 
tions obtained ratios measured between 4 and 50 GeV with a series of balloon flights. The 
CAPRICE data do not show any flattening of the ratio with increasing energy (ratio ~ ICR 3 
around 40 GeV) as expected from the secondary production model, in contrast with the last 
HEAT data point, which sets a limit of ~ 2 • 10” 4 above 20 GeV. This underlines the limited 
statistics available up to now in this kind of experiments. 

The uncertainties of theoretical models are large below the well defined flux maximum of 
2 GeV due to the complexity of production and propagation at low energy. On the contrary, 
the secondary high-energy flux is predicted with good confidence, all estimates being consistent 
with each other. Above a few tens of GeV, the antiproton production becomes quite negligible, 
the flux falling by 3 orders of magnitude below the maximum for antiproton energies around 
40 GeV. Any experimental hint of antiprotons in these high energy regions would therefore be 
of prime importance. 

First upper limits [12] on the p/p ratio around 1 TeV were presented by the L3+C and the 
TIBET-ASy - collaborations [13] at different conferences. The Tibet array has a worse angular 
resolution compared to the L3+C experiment and observes a smaller deviation of the Moon 
shadow due to its sensitivity to higher primary energies. 

In reference [14] an upper limit on the p/p ratio around 1 TeV is calculated from different 
measured /i + //i~ - ratios at ground level with large uncertainties. This indirect determination 
of a limit is based on cascade calculations which depend on the assumed primary composition 
and the hadronic interaction cross-sections at high energies [15]. The method has therefore 
large systematic uncertainties. L3+C has recently measured precisely the muon momentum 
spectrum, as well as the charge ratio and the angular dependence [16]. Based on to-day’s 
knowledge of the parameters entering the calculation, an estimate of an upper limit on the 
contribution of antiprotons to the primary flux could not compete with the one presented in 
this paper. 

There are at least 3 models of exotic sources able to produce high energy antiprotons: 
primordial black hole (PBH) evaporation [17], dark-matter neutralino annihilation [18,19] and 
high-energy antiprotons from extragalaetic sources [14,20]. In some scenarios, the last two mod¬ 
els can provide a p/p ratio increase up to the 10 % level in the energy range under investigation. 
This shows the importance of measuring the antiproton flux at high energy. The antiproton 
energy-spectrum is of course expected to be different for each type of sources. Putting a limit 
to the number of antiprotons will, in any case, constrain some of the parameters of the models. 


Moon shadow experiments EAS arrays were the first detectors to look for a Moon-shadow 
effect and the first observation was reported by the CYGNUS collaboration [2], with a 4.9 
standard deviations (s.d.) significance. The CASA group [3], with a larger array, obtained 
about the same result. In absence of any dominant point-source, the Moon-shadow experiment 
provides a unique possibility to make a direct measurement of the angular resolution of the 
detectors and to verify the pointing accuracy. 
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Large underground detectors also have the potential to observe the Moon-shadow. Back¬ 
tracking muons to the surface depends on the correct estimation of multiple Coulomb scattering 
in the large rock overburden, a process known to be essentially non-Gaussian. Results have 
been presented by the MACRO [5, 6 ] and SOUDAN [7] collaborations. These detectors have 
muon energy-thresholds of several TeV. As a consequence, statistics are low even if the data 
were accumulated for nearly 10 years. The effect of the geomagnetic field is hardly visible 
being less important than the observed shift from the origin due to pointing uncertainties. To 
make the influence of the deflection due to the magnetic field significant, the only possibility 
is to lower the detection threshold of the primary particles. Locating EAS arrays at very high 
altitudes is a solution. The TIBET air-shower experiment [4] has been operated since 1990 at 
4300 m above sea level ( a.s.l .) and has provided the first unambiguous effect of the geomagnetic 
field on the Moon shadow. 

Imaging-Cherenkov detectors have also been proposed for the observation of the Moon 
shadow [ 8 ]. The search of 7 sources with this technique was a success, in particular with the 
observation of the Crab nebula and a handful of other point sources. The application to the 
Moon-shadow measurement is more difficult as moonlight prohibits the use of visible photons 
and no Moon shadow was observed with this technique [21], 

Another promising technique uses a large volume of water as the detection medium. Photo¬ 
multiplier tubes detect the Cherenkov radiation produced in the water by relativistic charged 
particles or photons produced in the primary shower. The MILAGRO collaboration [22, 23] 
built a first prototype, MILAGRITO, running from February 1997 to May 1998, then a full 
detector, MILAGRO, starting its operation in February 1999. The goal is to be sensitive to 
primary cosmic rays down to 1 TeV or less, as imaging-Cherenkov detectors, while maintaining 
an all-sky acceptance and a high-duty cycle like EAS arrays. Preliminary results have shown 
that the Moon shadow effect is observed with a significance above 20 s.d. [24], No result of an 
antiproton search is yet available. 

2 The L3+C detector 

The L3+C detector is part of the L3 apparatus [25], one of the four particle detectors installed 
on the LEP Collider. It is located underneath the French-Swiss border at CERN, at 450 m a.s.l. 
under 30 m of sedimentary rocks called molasse (density ~ 7.2 x 10 3 g/cm 2 ). It mainly makes 
use of the muon chamber system which was designed to make a very precise measurement of 
muons produced in e + e _ collisions. The muon spectrometer consists of two octagonally shaped 
rings, each with eight “octants”, installed in a 12 m diameter solenoidal magnet which provides 
a uniform field of 0.5 Tesla along the e + e _ beam direction. Each octant contains precision drift 
chambers organised in three layers to measure the projection of the muon trajectory onto the 
plane orthogonal to the magnetic field, and layers of drift cells to measure the projection along 
the magnetic field direction. Other parts of the L3 detector are not used by L3+C. To fulfil the 
specific features of the cosmic ray experiment and to make the running of both L3 and L3+C 
completely independent from each other, several systems are specifically added to the L3 setup: 

• On top of the magnet, 202 m 2 of plastic scintillators are installed to determine the muon 
arrival time. 

• A new trigger and data-acquisition system is built to decouple the L3 and L3+C operation. 
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• A precise timing system is devised. It is based on an external GPS module and it includes 
also 1Hz and 10 MHz clocks. 

• An air-shower scintillator array is installed on the roof of the surface building to estimate 
the shower size associated to a detected muon. Its data are not used in the present 
analysis. 

The geometrical acceptance of the detector amounts to about 200 m 2 sr and the muon momen¬ 
tum threshold set by the overburden is 15 GeV. The detector was operational in May 1999 and 
a total of 1.2 x 10 10 muon triggers were collected up to November 2000, corresponding to an 
effective live-time of 312 days. Both momentum resolution and detection efficiency are checked 
using muons from Z decays that L3+C could also detect when the accelerator was running at 
a centre-of-mass energy equal to the Z-boson mass. Studies are extended using muons going 
through two separate octants, both giving a nearly independent measurement of the particle 
momentum. The momentum resolution is found to be 4.6% at 45 GeV and 7.4% at 100 GeV. 
The muon momentum threshold can be adjusted off-line to optimise the results. A detailed 
description of the L3+C detector and its performances is given in References [16,26]. 


3 Experimental considerations and Monte Carlo simula¬ 
tions 

A specific description [27] of the Earth-Moon spectrometer system and of the cosmic ray shad¬ 
owing effect is implemented in the L3+C simulation with two aims: 

• to take into account as accurately as possible the different detector components, the 
Earth magnetic held, the cosmic ray showering in the atmosphere, the multiple scattering 
of muons in the molasse, and the reconstruction, 

• to understand the relative importance of the various parameters contributing to the mea¬ 
surement of the p/p ratio. 

The detector properties are described with a Monte Carlo based on the GEANT program [28] 
which takes into account the effect of energy loss, multiple Coulomb scattering and showering 
in the detector. The basic version of the model is identical to the L3 simulation package, but 
specific features required by the L3+C setup, are taken into account, such as the additional 
scintillators on top of the magnet, or the magnetic field in the coil and yoke of the magnet. For 
the measurement of cosmic rays originating from the atmosphere, the overburden above the 
detector must also be included. The energy loss of muons and the smearing of their angular 
direction is an important issue. The whole surrounding of the L3 detector, consisting mainly of 
molasse, is introduced, including the access shafts to the experimental cavern and the shielding 
structures. All main physics processes related to the muon propagation through matter, such 
as multiple scattering, secondary-particle production including h-rays, pair production, energy 
loss and decay are fully simulated. 

Special attention is put on the simulation of the muon chambers, by including all inefficien¬ 
cies due to less efficient and dead cells in the muon detector. The simulated Monte Carlo events 
are reconstructed and backtracked to the ground level in the same way as the data events. 

The simulation of the detector is based on the generation of muons. Instead of performing 
a full simulation of the air-shower cascade generated by the primary cosmic radiation in the 


5 



atmosphere, single and double muons are generated above the detector, according to the known 
angular and energy distributions obtained by a full air shower simulation using the CORSIKA 
package [29]. The interactions, decays, annihilations and secondary-particle production in 
the air are fully simulated, according to the current experimental knowledge and to various 
theoretical models. 

3.1 Angular resolution 

One of the key issues of the Moon-shadow measurement is the experimental angular resolution. 
The Moon subtends a radius of 0.27° and the angular resolution has to match this constraint. 
The following factors are taken into consideration: 

• the muon direction with respect to the primary nucleon direction, 

• the multiple scattering in the molasse above the detector, 

• the intrinsic angular resolution due to the muon chamber resolution, the alignment and 
the reconstruction precision. 

As a result, the angular resolution is a complicated function, depending not only on the 
muon momentum and on the amount of matter on the particle trajectory, but also on the 
variables used in the event selection. A good opportunity to characterise the angular precision 
and check the simulation is given by the study of the space-angle distribution of two-track events 
in the detector, called “di-muon events” in the following. Muons coming from the decays of 
mesons originating from the early stages of the shower development are produced very high 
in the atmosphere and move along nearly parallel paths. The angular separation of the two 
muons is therefore a good measure of the smearing introduced in their direction by all the 
effects mentioned above. The results of the L3+C simulation is compared with the obtained 
di-muon data. 

3.1.1 Di-muon analysis 

Di-muon events are selected with cuts requiring a minimum quality of the two tracks. The main 
goal of the selection is to remove fake di-muon events i. e. events with single muon split into two 
different tracks because of reconstruction problems. For this purpose, a minimum separation 
between both tracks is required. Events are further classified into “double-double” “double¬ 
single” and “single-single” according to the number of subtracks for each of the reconstructed 
track, a subtrack being defined for each octant crossed. A muon momentum threshold, dehned 
at the ground surface level, is also imposed on both muons. Except for high muon momenta, 
large statistics are available and results are mainly dominated by systematic uncertainties. 
Events simulated with the same sets of cuts are compared with experimental results. Figure 1 
is an example of such a comparison for all events with muon momenta between 50 and 60 GeV. 

In the following, we define the di-muon angular resolution as cr 2/i = HWHM/1.17 with 
HWHM being the half-width at half-maximum of the distribution peak. The observed and 
expected values of cr 2/J are shown in Figure 2 as a function of the muon energy for the whole 
di-muon sample. 

Another check of the detector simulation follows from the study of the angular resolution 
versus the amount of matter crossed by the muons before reaching the detector. The multiple 
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scattering is the main factor in the contributions to the angular resolution from different com¬ 
ponents. A large range of matter thickness is available by selecting events from the access-shaft 
direction (minimum energy loss) or from large zenith angles (maximum energy loss). 

Experimental results are compared with the detector simulation results in Figure 3. In all 
cases, except at large angles, the data is in a rather good agreement with the results from the 
simulation. 

3.1.2 Angular resolution in the Moon shadow analysis 

The present analysis uses single-muon events whose angular resolution cannot be deduced from 
di-muon data. In addition, as the distribution of both type of events inside the air shower are 
quite different, they do not have the same acceptance. The good agreement observed between 
data and Monte Carlo in the di-muon analysis gives confidence in the use of the simulation to 
extract the angular resolution in two ways: 

• The angular resolution can be directly extracted from the event deficit in the Moon 
direction. As already mentioned this is a unique opportunity for an angular resolution 
measurement and the results will be compared with the simulation. 

• Alternatively, the angular resolution can be fully constrained with the help of the sim¬ 
ulation of single muon events. The point-spread function for each momentum bin is 
determined using the same track selection and the same acceptance as for the selected 
Moon events. The corresponding angular distributions can then be used in the shadow 
simulation as smearing sources for the angular resolution. Systematic uncertainties are 
at the level of 5 % or less, better than the angular resolution which can be extracted from 
the di-muon data. These systematic uncertainties are estimated from the comparison of 
data and simulation in the di-muon results and from the studies of different production 
models in the simulation of air showers with the CORSIKA program. 

3.2 Earth magnetic-field and calculations of the deflection 

Two models are investigated to describe the geomagnetic held, a simple dipole model and the 
International Geomagnetic Reference Field model (IGRF) [30]. In the latter, the geomagnetic 
held is commonly expressed as the gradient of a scalar potential which can be expanded in 
terms of spherical harmonics. The IGRF consists of a series of values of the coefficients in 
the expansion based on direct measurements of the geomagnetic held. In fact, the hrst terms 
of the expansion can be identified with the held produced by a dipole located at the centre 
of the Earth. The contributions of the other terms can be considered as perturbations of the 
main dipole held. Quantitative differences between results of both models have been studied. 
The main difference is a shift of the southern magnetic pole. As a consequence, differences at 
the 10 % level are observed at the L3+C location concerning the held intensity. Differences on 
the magnetic held direction can reach 5° in the part of the sky where the Moon is visible by 
the detector. As a result, calculated dehections can differ at the same level. The amount of 
deflection is overestimated by the dipole model in the largest part of the sky and reaches 10 % 
for a 1 TeV proton. Consequently only the IGRF model is used in the following. 

Two coordinate systems are used. The hrst is based on the local horizon. The zenith angle 
9z and the azimuth angle ccz are determined from the precise position of the detector. The 
second is an equatorial system using the Earth rotation axis as basis for the two coordinates: 
declination, <5, and right ascension, RA. 
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The deflections of the particle trajectories in the magnetic field is described in an additional 
coordinate system. During its way from the Moon to the Earth, a particle of charge Z and 
momentum p is subject to the Lorenz force in the field B and the deflection A 9 is linked to 
the particle path /, as: 


A#(mrad) 


°' 3 ' p(TeV) 


B (Tesla) x d l (m) 


( 1 ) 


Depending on their incident direction at the top of the atmosphere, charged particles traverse 
different field regions. Thus, for a given particle, the angular deflection is a function of the 
incident direction, the charge and the momentum of the particle. This can be used to establish 
a deflection map that gives, for a given momentum, the amount of deflection and its direction. 
Figure 4 shows a Moon transit above the sky as seen by the L3+C detector. Each point 
corresponds to one direction in the sky as computed from the zenith and azimuth angles. 
During a Moon transit in the sky, the direction of deflection strongly depends on the Moon 
position but the dependence on the momentum is rather small. This leads to the definition 
of a coordinate system defined for each Moon position in the sky, with coordinates 6 *h and 
9y respectively parallel and orthogonal to the direction computed for a particle with a given 
primary momentum (here a 1 TeV proton). The indexes H and V stand here respectively 
for ‘horizontal’ (parallel deflection) and ‘vertical’ (perpendicular to deflection). In this way, 
magnetic deviations will shift the Moon shadow image along the parallel direction and the 
shape in the other direction will mainly depend on the angular resolution. 


3.3 Primary cosmic-ray composition 

Energy spectra for various elements, up to a few hundreds TeV for protons and a few tens 
TeV per nucleon (TeV/N) for heavy elements, have been measured with the use of balloon and 
satellite experiments. The proton spectra obtained by different experiments are in reasonable 
agreement. Results are considerably scattered for other elements, a consequence of limited 
statistics and uncertainties in the energy calibration. 

A compilation of available data [31] proposes the following fit for the flux of particles: 

$ = $ 0 E~^ A ( 2 ) 

with E the energy per nucleus in TeV. The power index for He 4 is smaller than the one for 
protons, and therefore the contribution of a particles increases at high energy. However, recent 
results from the RLINJOB [32], AMS [33] and BESS [34] collaborations seem to invalidate such 
behaviour, with a common power index ~ 2.8 for proton and helium spectra. 

The composition of primary cosmic rays plays an important role for the Moon shadow. It 
acts remarkably differently for experiments using EAS arrays, Cherenkov or /j detectors. For 
the first two methods, the measured signal is proportional to the total energy E of the primary. 
The third method, characterised by the /j momentum threshold, is sensitive to the primary 
energy per nucleon E n = E/A. 

Muons with energy E M are produced by nucleons of a minimum energy Eq with Eq ~ E lt 
and thus by nuclei with energy E > AE 0 . In a nucleus, all A nucleons may contribute to the 
interaction. Let us call <j(E A N ) the cross section for the production of a muon with energy 
Ep by a nucleon with energy E^. If the spectrum has an index 'Ja and the corresponding flux 



( 3 ) 


is (pA at 1 TeV, then the number of muons with energy E [L produced by these nuclei is: 

roc 

N a (EP) = A aiE^E^E-^dE. 

JAEo 

If the energy spectrum of all particles follows a power law with the same index 7 and the 
probability to yield a muon with energy E M does not depend on the energy per nucleon above 
threshold, the contribution of a nucleus with A nucleons relative to the proton contribution is : 


N a {E„) 

ME,) 


oc r A A 7+2 , 


(4) 


where r a = 77 is the relative abundance of a nucleus compared to the proton abundance at 
1 TeV. Using data from Reference [31], it is found that muons originate at 75.8 % from protons, 
17.3% from helium nuclei and 6.9% from heavier nuclei. 

According to equation ( 1 ), for a mean primary energy Z mean , the mean deflection angle is 
proportional to the ratio —. For experiments sensitive to the total energy, like EAS arrays, 
the mean primary energy does not depend on A. Thus 


(Ad) eas oc Z (5) 

The position of the Moon shadow depends only on Z, therefore EAS-array experiments 
expect separate shadows for protons, helium and heavier nuclei. 

For a muon experiment, the mean primary energy is proportional to AE 0 . Thus the mean 
deflection angle is: 

(A0) m oc ^ (6) 

Muon experiments are sensitive to the ratio which is equal to 1 for protons, and from 0.5 
to 0.4 for heavier nuclei. All shadows from helium and heavier nuclei are almost at the same 
place. Therefore, in the following, the primary flux for the observed muons will be considered 
to be 75 % protons and 25 % helium nuclei. 


3.4 Primary cosmic-ray energy spectrum 

There is an energy window for the observation of a magnetic-field effect on the Moon shadow. 
High primary energies (> 10 TeV) allow ground level observations with large detectors and 
relatively good statistics. However the magnetic deflection will be small compared to the 
angular resolution, making the effect negligible or, at most, appearing just as a small correction. 
Low primary energies (~ 100 GeV) are difficult to observe from the ground and large deflections 
due to the Earth magnetic held will dilute the shadow image and severely limit the sensitivity. 
L3+C has a good sensitivity to muons from low energy primaries. Moreover, the muon energies 
are measured with good precision. The shadow effect can be observed using different ranges of 
muon energies, thus selecting samples of different primary energy spectra. For each observed 
muon energy E a corresponding primary energy E with E > E [L is obtained using the shower 
generation with CORSIKA and the tracking of muons with the detector simulation. Figure 
5 shows the expected proton and helium spectra associated with a detected muon with E = 
100 GeV. The maximum of the primary energy distribution is around 1 TeV for protons and 
4 TeV for helium nuclei. 
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3.5 Moon-shadow simulation 


The simulation program tries to reproduce as closely as possible the conditions and the param¬ 
eters which significantly influence the observation of the Moon shadow. However to make the 
simulation more efficient, the particles are assumed to be coming from the Moon surface: a 
positive signal is simulated instead of a deficit. Also, for the same reason, particles are followed 
backward instead of forward. They are originating from the detector and tracked through the 
Earth magnetic field up to the Moon. The momentum and identity of the primary particle are 
extracted from the distributions described above. For the angular smearing two methods are 
applied: either an arbitrary effective angular resolution is used, or the momentum dependent 
angular information coming from the experimental detector simulation. Events are generated 
during a time span reproducing the experiment running time and the detector acceptance. 
Some examples of simulation results are shown in Figure 6. 

The shape of 6 h becomes more symmetrical as the angular resolution is worsening. In the 
orthogonal direction, the smearing contribution is mainly coming from the angular resolution 
only. The expected signals for primary protons and helium nuclei are largely overlapping and 
no extraction of the helium contribution can be expected from the data. Instead of the p/p 
ratio (rp/ p ) mentioned before, the analysis will try to get the “p content” as seen by L3+C, 
r = 4>p/(j} matter, with </> mat t er the flux responsible of the observed deficit and assuming no anti- 
Helium contribution. The p/p ratio itself can then be deduced from this result and the estimated 
proportion (75 %) of the deficit due to the proton flux relative to the the total matter flux, as 
discussed in section 3.3. 

Figure 7 shows the simulated Moon shadow as it appears for E / > 100 GeV in the three 
coordinate systems: local, equatorial and deflection. Both the offset and the elongation due 
to the magnetic field are more visible in the deflection system. The search for a possible 
“anti-shadow” due to antiprotons is therefore performed in this last system. 

There are two ways to use the simulation results in this analysis, either directly, or through 
a parametrisation. In the first case, simulated distributions represent the expectation values to 
be compared to the experimental data. In the second case, the adjustment of the parametric 
function to the data allows the extraction of the parameters. This is used to give a result on 
the observed effective angular resolution. The simulated shadow is parametrised as the product 
of two functions, one for each direction: 


f(x,y) = fn(x) x f v (y), 


(7) 


/v can be described as the projection on the vertical axis of the two dimensional convolution of 
a disk (the Moon) with a Gaussian distribution, corresponding to an effective angular resolution 
for the muon momentum range considered: 


fv(y) 
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( 8 ) 


For /h, no analytical description exists. However, it is found that a sum of two Landau 
distributions approximates reasonably well the shape in the direction parallel to the deflection: 


fn(x) = (1 - c) • Li(x, Q'i, h ) + c • L 2 (x, a 2 , b 2 ), (9) 

with the parameters a and b acting respectively on the width and the position of the maximum 
of the distribution. For each selected value of the angular resolution, a number of events corre¬ 
sponding to 100 times the number of expected deficit events in the experiment is simulated. A 
good agreement between the simulated shape of the deficit and the parametrisation is observed. 
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4 The Moon data-samples 

To define the muon direction in the sky, both local and terrestrial based coordinate systems 
are used. The Moon position is computed using the “SLALIB” library subroutines [35]. The 
error on the position calculation is estimated to be smaller than 0.01°, much smaller than the 
angular resolution or the Moon radius. The Moon is the nearest and the only astronomical 
object for which the position is significantly dependent on the observation location on Earth. 
This parallax effect is taken into account. It amounts to a few tenths of a degree. 

Another effect is the change of the apparent size of the Moon as seen from the Earth due to 
the variation of the centre-to-centre distance from the Earth to the Moon between perigee and 
apogee. Calculations show that the Moon angular radius, as seen from the Earth, fluctuates 
between 0.25° and 0.28° during the whole data taking period. 

The Moon declination is continuously changing inside a range of around ±20°. In the local 
sky, the Moon follows a trajectory reaching a minimum zenith angle of 25° for the experimental 
running period. 

As the detector cannot be triggered above a certain zenith angle value, the Moon is only 
available for certain periods of time, each called a “cycle” in the following. Figure 8 shows the 
Moon acceptance for the selected events for the two running years. Data were accumulated 
for Eve cycles in 1999 (73 transits) and nine cycles in 2000 (142 transits). The Moon was 
available for a total of 1557.5 hours. The corresponding data-acquisition live-time is 1188.7 
hours (76.3%). 

Data selection and monitoring A first selection isolates events coming from the direction 
of the Moon. For each event, the Moon position is computed in local coordinates and the space 
angle 6 with the muon track direction is deduced. The track-reconstruction program requests 
at least one “triplet” (hits from 3 chambers in one of the octants) and one scintillator hit. Two 
hits in one octant (called a “doublet”) are also allowed if a momentum mesurement is possible, 
but two doublets are rejected. This constraint leads to the reconstruction of only one third of 
the total number of events collected. The following cuts are applied at this level: 

• Only one muon is reconstructed in each event. 

• The angle with the Moon direction is less than 5°. 

• The muon momentum is more than 50 GeV. This cut is motivated by the fact that low 
energy muons have little correlation with the direction of the primary cosmic ray. 

• The Moon zenith-angle is less than 60°. Above this value reconstruction becomes more 
difficult and the trigger efficiency is low. 

• The backtracking of the measured muon track in the detector up to the surface is suc¬ 
cessful. 

• Events with timing uncertainties, amounting to 0.2% of the total, are rejected as correct 
time information is needed to compute precisely the Moon position. 

When these cuts are applied, a total of 6.71 x 10 5 events are selected, out of which 2.11 x 10 5 
in 1999 and 4.60 x 10 5 in 2000. 

To monitor the data, some variables are carefully tested as a function of time. Among them 
are the number of selected events, the number of high energy events, the /i + /charge ratio, 
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the proportion of muons including 2 subtracks and the proportion of high quality subtracks 
in events. Di-muon events are also used. They are taken from directions along the Moon 
trajectory using the same runs as in the Moon data. A total of more than 50000 events are 
collected in this way. The event selection requires that at least one of the muons has 2 subtracks 
and that the minimum muon momentum is 50 GeV. The di-muon space angle distribution is 
monitored. 

Stability with time is an indication that no major problem perturbed the collection of data 
during all the running of the experiment. No major instability is observed for the whole running 
period. 

Background determination Searching for a point source (or shadow) entails the counting 
of the number of events in an angular bin containing the possible signal (signal bin) and 
comparing it to the number of background events expected in this bin. The background is a 
function B(9 z ,a z ,t) of the zenith angle 9 Z , the azimuth angle a Zl and the time t. 

Apart from the 9 Z dependence due to the changing thickness of the atmosphere, there is 
a strong spatial dependence of the reconstruction efficiency clue to the arrangement of the 
muon chambers in octants and to the constraints imposed on the track reconstruction. The 
background is therefore evaluated by counting events in regions that were exposed for the same 
amounts of time to the same directions of the sky as the signal bin. 

The global rate is changing with time, due to modifications in the detector hardware or 
in the detector environment (noise dependence, local atmospheric temperature and pressure 
dependence). A correction has therefore to be applied. In general time and spatial angular 
dependence are independent and the spatial acceptance is nearly constant. If not, one has to 
consider sufficiently small time slices so that this is valid. 

The background is determined by measuring the number of events due to “fake Moons”. 
These are bins which cross a given region in the sky either earlier or later than the signal 
bin in different runs. Averaging background samples on both sides of the signal bin removes 
effects of changes in the event rate which are linear in time. Four “fake Moons” are used, one 
hour and two hours before and after the real Moon position. When considering both running 
years, ten samples are available for the background evaluation (the two signal samples and eight 
background samples). 

The projections on the azimuth and on the zenith axes are shown in Figure 9a and Figure 
9b. In these examples, the merged data from all samples are used. The distributions are fitted 
by a straight line. The fitted slope parameters from the individual samples are shown in Figures 
9c and 9d for the azimuth and zenith directions respectively. There is only a slight positive 
variation of the rate as function of the azimuth angle and all the samples give statistically 
compatible results. The variation with the zenith angle is more important. The rate is de¬ 
creasing for large zenith angles, a consequence of the acceptance. Moreover samplc-to-samplc 
fluctuations are larger. 

Each sample is fitted by a plane in the azimuth-zenith coordinate system. The expected 
background density at the origin (real or “fake Moon” position) is computed, expressed by the 
number of events in the Moon solid angle, Dmooii, for the given live-time of the whole experiment. 
Fluctuations between samples are at the level of a few percent, much greater than the statistical 
uncertainties. Therefore systematic uncertainties clue to live-time errors and acceptance or rate 
changes with time dominate. Following the hypothesis mentioned above that time and spatial 
angular dependence are independent, a further normalisation correction-factor based on the 
overall number of events for each sample can be applied. This is obtained from the total 
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number of events inside an annulus around the nominal Moon position (3° < 6 < 5°). After 
this normalisation, differences are at the 0.1% level. When applying the above procedure, 
there is no significant difference in the evaluation of the background using different samples. 
The systematic uncertainties in the knowledge of the background rate at the Moon position 
are negligible compared to the statistical uncertainties on the signal. The result averaged 
for the two years is 542 ± 0.6 events/ffMoon — 2366 ± 3 events/deg 2 for E / > 100 GeV and 
677 ± 1.0 events/flMoon — 2956 ± 4 events/deg 2 for 65 GeV < E tl < 100 GeV. 


5 Event-deficit analysis 

In the local coordinate system, evidence for a cosmic ray deficit introduced by the Moon is 
observed using a single angular variable. The number of events shows a linear increase with 
increasing solid angle. Therefore the density of events is generally considered. In absence of 
any signal, plots must show a flat distribution. The plots of Figure 10 show the results for 
E, > 100 GeV with a “fake Moon” shifted 1 hour behind its real position along its trajectory 
(Figure 10a) and with the Moon at its nominal position (Figure 10b). In this last plot, a clear 
deficit of events in the first few bins is observed. This is attributed to the shadowing effect of 
the Moon. 

Uncertainties in pointing and the influence of the geomagnetic held are also contributing 
to the shape of the deficit. The extraction of the experimental angular resolution and the 
measurement of the p content in primary cosmic rays are not possible on this one dimensional 
distribution. In the following, a maximum likelihood method is used to disentangle the various 
factors. An interesting property of the “deflection coordinate system”, which is based on the 
deflection induced by the Earth magnetic held, is to concentrate the Moon shadow deheit along 
one axis, thus optimising the signal density. Therefore the analysis will be performed in this 
coordinate system with the two projection angles and 0y. Studies to investigate the effect 
of the muon momentum threshold on the deheit lead to the definition of two samples, a “high 
energy (HE)” sample for E / > 100 GeV and a “low energy (LE)” sample for 65 GeV < E^ < 
100 GeV. 

Figure 11 shows “shadow” plots concerning the data for both samples. Smoothing techniques 
are used. The way the background is computed is described later. For the analysis, “raw” 
spectra are used. A binning of 0.1° is chosen in each direction. 

The shape of the shadow is more elongated in the case of the LE sample and its position is 
shifted further. The position and the shape of the Moon-related deheit mainly depend on the 
magnetic deflection undergone by the primary particle associated with each muon of the sample 
and on the effective angular resolution a. Both effects are taken into account in the simulation. 
The effective angular resolution includes the muon production-angle inside the air shower, the 
multiple scattering in the molasse above the detector and the intrinsic angular-resolution due 
to muon-chamber resolution, alignment and reconstruction. The smearing due to the multiple 
scattering is the main component of a. A parametrisation of the simulation output for each 
value of a and for each sample was described in 3.5. It is used in the maximum likelihood fit 
to allow the extraction of a as a free parameter. 

Analysis procedure The probability to have n l:] events in bin (i, j) of the 0y — Oh distri¬ 
bution, when g(xi,i/j) events are expected is discribed by Poisson statistics. The logarithm of 
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the likelihood function can be written as: 

N N 

in £ = ^2 [nijHgi^yj)) ~g{xi,yj) - ln KtO]- 

i,j =1 i,j =1 

For practical reasons, the likelihood function is normalised and a minimum A In £ m is searched 
as a function of the parameter values. The two-dimensional distributions are the result of the 
combination of 3 different components: 

• a smooth background, which can be fitted with a plane, 

• the proton and helium deficit introduced by the Moon shadow, 

• the antiproton deficit if any. 

The simulation provides a description of the proton deficit. A similar description is used 
for the antiproton deficit. However, in the deflection coordinate system, the shadow position is 
inverted with respect to the Moon centre. Also the shape of the antiproton shadow will differ 
from that of the protons, due to a possible different power index 7 of the energy spectrum. 
Thus the most general description of the data is: 

Amiss 

g(x,y) = u x x + u y y + u z --^-[0.75 fi(x - x Q ,y - y 0 , a) +0.25 f 2 (x - x 0 ,y - yo,<?) + 

'-v-' 1 + r "- v -' "- v -' 

background p deficit He deficit 

+ rf 3 (x 0iV -x,yo,p-y,(Tp)], (11) 

S --V-' 

p deficit 

where / 1 , f 2 and / 3 are the shadow functions defined by equation (7), respectively for protons, 
helium nuclei and antiprotons. The parameters that can vary for the maximum likelihood fit 
are the angular resolutions a, the positions of the deficits (a; 0 , Vo), th e number of missing events 
Ami^ and r, the p content. The parameters u x , u y , u z describe the background. 

The influence of different power indexes 7 (—1.7 to —3.7) of the antiproton energy spectrum 
on the final limit of the antiproton to proton ratio presented in this paper is studied, as well 
as a re-analysis of the data assuming also different values of 7 . Variations of less than 20% 
with respect to a simplified model, where the power index is assumed identical for protons and 
antiprotons are obtained. The interpretation of this conclusion may be explained by the fact 
that for a steep antiproton spectrum the deflection of most antiprotons is relatively large, but 
the muon statistics small, due to the lower average primary energy. In the case of a flatter 
spectrum the deviation is small, but the muon statistics larger. The two effects cancel each 
other, providing a negligible influence on the limit of the p/p ratio. For simplicity we present 
only the result of the analysis with the assumption of equal power indexes, reducing thus the 
number of free fit parameters to 8 (/ 3 = / 1 , £ 0 ,p = — x 0 , y 0t p, = — y 0 , 07 = cr). 

Instead of trying to extract directly the 8 parameters with the simplhed equation (11), the 
analysis proceeds in several steps. 

Background estimation A first determination of the background parameters is performed 
using the ring-data defined as 3° < 9 < 5°, where 9 is the angle between the muon and the 
nominal position of the Moon. The cut excludes the cells in the proton and antiproton deficit 
regions. The event density at the nominal Moon position is known at the 0.3 % level. This 
uncertainty corresponds to the statistics used for its determination. 

No significant changes in the parameter values are found when equation (11) is applied to 
the whole angular range and all the parameters are considered free. 
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Table 1: Results obtained in the fit of the matter deficit. 



HE 

HE 

LE 

LE 

parameter 

measured 

expected 

measured 

expected 

x 0 

0.33 ± 0.08 

0.26 

0.53 ± 0.13 

0.48 

Vo 

0.05 ± 0.05 

0.0 

- 0.10 ± 0.08 

0.0 

FWHM 

L07 -0.04 

1.03 

1.80 ± 0.15 

1.87 

N miss 

57511 ? 

546 ± 5 

5361^3 

683 ± 6 

Significance 

8.3 s.d. 

8.0 s.d. 

5.5 s.d. 

5.8 s.d. 


The Moon shadow analysis In this step no antiprotons are supposed in the primary flux. 
The effective angular resolution of the detector, together with the pointing uncertainties, are 
obtained from the observation of the matter deficit only. Results of the maximum likelihood 
fit for the two samples are given in Table 1. As an example, two-dimensional 68% and 90% 
confidence level contour curves for the parameters N m ; ss and a are shown in Figure 12a for the 
case of the HE sample. 

Pointing errors are given by horizontal or vertical offsets between data and simulation in 
the determination of the deficit position xq and z/q. Both values are small (< 0.1°). 

Values related to the absolute position, x 0 , and to the extension of the deficit (FWHM, full 
width at half maximum) in the horizontal direction $h show clearly a momentum dependence. 
In the 6\ direction, no shift is observed and the width is mainly the result of the effect of a. 

Values of V miss are extracted from the fit. However V miss can be also directly deduced from 
<F, the flux measurement around the Moon direction, V miss = x Ti ive x HMoon where Ti ive is the 
live-time corresponding to the Moon observation. The main contribution to the uncertainty 
is related to the Moon solid angle. At any time the precision of the calculations of the Moon 
radius Rm , from the SLALIB subroutines [35], is estimated to be 0.4 % and the time distribution 
of the Moon events is very well known. The uncertainty on the solid angle, proportional to 
R“l j is thus 0.8%. The flux is known at the level of 0.3% and the uncertainty on the live-time 
is still smaller. In total, the precision on the most probable value of N m j s s is estimated to be 
around 1%. The expected values of V miss are shown in Table 1. This knowledge of V miss is 
introduced in the likelihood function and allows an improvement in the determination of the 
remaining parameter a. N m - lss is constrained to its most probable value, with a 1 % uncertainty. 
The results are illustrated in Figure 12b for the LE sample with the plot of Aln£ versus a. 
The experimental results are a = (0.22 ±0.04)° for the HE sample and a = (0.28 ^ 05 )° f° r the 
LE sample. These numbers refer to an effective angular resolution valid for the set of selected 
events, integrated over the momentum distribution of the data and the directional range of the 
Moon events. For the HE sample, Figure 13 shows a comparison of a projected band of data 
around the nominal position of the Moon with the fitted results corresponding to Table 1. 

As a cross check, the effective angular resolution is also obtained by a fit to the one- 
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dimensional deficit distribution. The values discussed above are confirmed, albeit with much 
larger uncertainties. 

The antiproton search To set a limit on a possible cosmic ray antiproton component, 
the number of missing events 7V miss is supposed to be shared between protons, helium and 
antiprotons. The total number of missing events has been constrained to the expected value. 
The antiproton deficit is described like the proton deficit with the corresponding parametrisation 
function symmetric to the proton one with respect to the Moon position. 

As mentioned in section 3.1.2, the value of a can be deduced from the simulation. Com¬ 
parison with data in di-muon events shows a good agreement. Contrary to the situation in the 
Moon-shadow experiment, high statistics is available and detailed investigations are possible. 
The obtained results are a = (0.24 ± 0.01)° for the HE sample and a = (0.38 ± 0.02)° for the 
LE sample. Uncertainties are better than those obtained above. This angular information is 
used and a maximum likelihood fit is performed using the p content as a free parameter. 

The HE and LE results are combined to give the final measurement. This is done simply 
by adding the likelihood logarithmic-functions of each range. The total significance of the 
deficit is 9.4 s.d. The uncertainty range is obtained by finding the points for which A In £ = 
Aln£ m + 0.5. The result is shown in Figure 14 and one finds r = </>p/0 ma tter = —0.07 ± 0.09. 
The result is below a physical boundary (the p content must be positive). An upper limit of 
0.08 with 90% confidence level is set using the unified approach [36]. With the assumed flux 
composition around 1 TeV of 75 % protons and 25 % heavier nuclei responsible of the observed 
deficit, this corresponds to a p/p ratio of rp/ p = 0.11. Figure 15 shows the L3+C result together 
with other published values. 


6 Conclusions 

The L3 detector has collected more than 10 10 triggers of cosmic ray muons during the years 
1999 and 2000 in parallel with high energy particle physics studies at the LEP accelerator at 
CERN. About 6.7 x 10 5 events, with a direction pointing to a 5.0° cone around the Moon, are 
used and a Moon-shadow effect in cosmic rays is observed. 

A two dimensional analysis confirms that the alignment of the detector is correct to better 
than 0.2° and that the size and the shape of the deficit are compatible with the expectations. 
Two sets of data corresponding to high ( E / > 100 GeV) and low energy muons (65 GeV < 
E m < 100 GeV) lead to values of the effective angular resolution respectively of (0.22 ± 0.04)° 
and (0.28 /[j'Jjj) 0 . These numbers include all effects due to the showering of the primary cosmic 
ray in the atmosphere, the multiple scattering in the molasse and the detector resolution. They 
describe correctly the observed event deficit. The observed significance of the Moon shadow 
effect is 9.4 s.d. A significant effect clue to the Earth magnetic field is observed. This is better 
seen in a coordinate system with axis respectively parallel and orthogonal to the deflection 
defined for each direction in the local sky. The offset and the extension of the shadow are clearly 
dependent on the muon momentum range considered. With the hypothesis that the presence 
of antiprotons in cosmic rays would lead to a symmetric shadow to the one due to protons, a 
measurement of the p content is extracted from the data and is found to be r = —0.07 ± 0.09. 
A 90% confidence level of 0.08 is set on r, corresponding to an antiproton over proton ratio of 

rp/p = 0.11. 
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Figure 1: Distribution of the square of the angle between muons, for data and Monte carlo. The 
distributions are normalized to unit area. 
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Figure 3: Di-muon angular resolution versus matter thickness. 






sin(0 z ) • sin(a z ) 

Figure 4: The Moon passing through the L3+C acceptance. The acceptance is determined by the particular 
structure of the drift chamber and scintillator assembly. The contour lines correspond to an observed cosmic-ray 
flux of respectively 75%, 50% and 25% of the maximum of the flux. A Moon transit is indicated with dots. For 
each dot, the geomagnetic-deflection direction and amplitude for a 1 TeV proton is indicated by an arrow. 
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Figure 5: Distributions of proton and helium primary energies for = 100 GeV. The normalisation between 
the distributions follows the chosen ratio: 75 % of protons and 25 % of helium. The vertical scale is in arbitrary 
units. 
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Figure 6: Examples of results of the simulation of the shape of the Moon shadow: a) $h for E M > 100 GeV 
and for several hypothesis for the angular resolution, b) 9y for E M > 100 GeV and for several hypothesis for the 
angular resolution, c) 0 h for three muon energy ranges, d) On for proton and helium primaries for E M > 100 GeV. 
All distributions are normalised to the same number of events. 
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Figure 7: Contour map of the simulated Moon shadow for E M > 100 GeV: a) in the local coordinate system, 
b) in the celestial coordinate system, c) in the deflection system. For each system a total of 6 • 10 4 events have 
been simulated. The solid contour line defines the domain containing more than 700 events per sky unit cell 
(of size 0.1° x 0.1°), and the dashed and dotted lines lower minimal numbers. The deflection system obviously 
records the largest amount of cells containing at least 700 events each. 
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Figure 8: Moon acceptance for a) year 1999 and b) year 2000. The plot shows the event arrival-time in 
days and hours. The running period includes 5 lunar cycles in 1999 and 9 in 2000. White bands during cycles 
correspond to periods in which either the detector or the data acquisition were not operational. 




































































Sample Sample 


Figure 9: Number of events versus a) the azimuth and b) the zenith direction, using the merged data from 
the 10 samples. The origin is taken at the real or “fake Moon” position, c) Slope parameters extracted for the 
10 considered samples in the azimuth direction and d) in the zenith direction. The first 2 bins are from signal 
and the remaining from background. 
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Figure 10: Angular distributions for events/Moon-solid-angle with > 100 GeV, (a) for a “fake Moon” 
shifted 1 hour behind its real position, (b) using the correct Moon position. Solid lines are the results of the 
simulation, including the angular resolution deduced from the study of di-muon events. 
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Figure 11: Results obtained in the deflection system for: a) the High Energy sample, b) the Low Energy 
sample. In both cases, smoothing techniques have been applied. A circle indicates the true position of the 
Moon. The vertical grey scale shows the significance in standard deviation units; negative values correspond to 
an event deficit. 
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Figure 12: a) Two-dimensional 68% and 90% C.L. contour curves for the parameters iV m i ss and er for the 
HE sample, b) Aln£ versus a for the LE sample. N m ; ss has been constrained. The dashed line for Aln£ = 0.5 
is used to determine the 68.3% central confidence interval. 
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Figure 13: Comparison between data and fitted results for the HE sample in a projection on a) the horizontal 
axis, b) the vertical axis. Differences in the projected bandwidths (a: 0.8 degrees, b: 3.0 degrees) explain the 
differences of the observed number of events. The dotted lines represent the fitted average muon event number 
outside of the Moon region, and the dashed lines the fit to the data in the Moon region. 
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Figure 14: Aln£ as a function of r (the p content) for the whole data. The dashed line Aln£ = 0.5 is used 
to determine the 68.3% central confidence interval. 
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Figure 15: Measurements of the ratio of the antiproton and proton fluxes versus the primary energy, including 
the L3+C limit around 1 TeV. (The dashed lines show the range of the theoretical expectations [9].) 
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